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Abstract 

In this work we discuss the natural appearance of the GeneraHzed Brackets in systems with 
non-involutive (equivalent to second class) constraints in the Hamilton-Jacobi formalism. We show 
how a consistent geometric interpretation of the integrability conditions leads to the reduction of 
degrees of freedom of these systems and, as consequence, naturally defines a dynamics in a reduced 
phase space. 
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1 Introduction 

Constraints in Lagrangian systems appear when the Hessian matrix built from the Lagrangian function 
is singular. Since the Hessian is the Jacobian of the Legendre transformation that takes velocities to 
conjugate momenta, the non existence of an inverse transformation leads to a phase space that is 
not isomorphic to the tangent bundle of the configuration space. In this phase space the symplectic 
structure is degenerate and a unique Hamiltonian cannot be defined, leading to problems in the usual 
construction of the Hamiltonian formalism. This problem was first treated by Dirac's [Tl|2l[3] canonical 
approach, in which the constraints are added to the canonical Hamiltonian function with Lagrange 
multipliers. The so called Dirac's algorithm is based on consistency conditions, that imposes the 
constraints as invariants over time evolution. 

Other alternative approaches have been formulated since, including the Hamilton-Jacobi (HJ) ap- 
proach, first developed in [3] and improved later in O El [7|, based on Caratheodory's equivalent 
Lagrangian method [8]. The Caratheodory's approach to HJ formalism has a very profound geomet- 
rical meaning: the necessary and sufficient condition for an Action to be minimized in a region of the 
configuration space is the existence of a family of surfaces everywhere orthogonal to a congruence of 
curves in this space. The family is a solution, in the case of regular systems, of a single HJ partial 
differential equation (HJPDE). If the system is singular with k constraints we will not have one, but 
a set of k + 1 HJPDE. The solution of the characteristics equations related to the HJPDE is the 
congruence of curves, called simply the characteristics, which forms a dynamical system with several 
independent variables. This geometrical interpretation is called in the literature [8] "the complete 
figure" of variational calculus. 

The integrability conditions are in the central role of the HJ formalism for singular systems. They 
should guarantee, supposing linear independence of the independent variables, the existence of a com- 
plete set of HJPDE and, by consequence, the integrability of the characteristics equations. The Probe- 
nius' integrability condition sets up that the constraints must form a complete system in involution 
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with the Poisson Brackets (PB). However, there are systems that do not obey this condition, presenting 
non-involutive constraints. To circumvent this problem some authors elaborate alternative methods 
[UllO], out of the scope of the HJ formalism. 

We will show that imposing the Frobenius' conditions in the form d{(j)} = 0, in which {i;^} represents 
the set of all constraints, we can deal with non-involutive constraints in a very natural way, allowing 
linear dependence of the former independent variables. The integrability conditions in this form also 
leads naturally to a non degenerate symplectic structure in a reduced phase space, represented by some 
Generalized Brackets (GB). The appearance of GB in the HJ formalism was already noticed in first 
order Lagrangians [11], but in this work we intend to make a more general formulation, to be valid for 
general, not only first order systems. 

2 The Complete Figure 

Let it be a Lagrangian function L(y,y,x,x,t), with n degrees of freedom, whose Hessian matrix is 
singular of rank m. This means that there are k = n — m constraints (/>^ = vr^ — dL/dx^ = and 
the system has its configuration space separated in two sub-spaces, namely the space of m variables y, 
relate to the invertible part of the Hessian, and the space of k variables x, related to the non invertible 
part. The equation 

vro + Pay" + vr^i^ -L = , 

where ttq = dtS, pa = daS and vr^ = d^S and the indexes are given by {a} = {l,...,m}, {z} = 
{1,...,A;} and {i} = {l,...,n}, is the sufficient condition for an extremum of the action J Ldt. 
Wherever the constraints are valid this equation becomes a HJPDE with the canonical Hamiltonian 
defined by 

Ho = Pay'' + TTzX^ -L , (1) 

which does not depend on the undetermined velocities x. Hence, the system must obey the set of 1 
HJPDE 

(/.„ = X„(5) = 7r, + i7„ = 0, {a} = {0,l,...,fc} , (2) 

in which = Xa^^i is a set of vector fields, Xa = ^^'^ Xa — d4>a/dpa- In our notation x^ is related 
to the time parameter and Hz = —dL/dx^. 
The characteristic equations are given by 

dy" = X«a!x" = K,M^^^" , (3) 

dPa = {Pa, 4>a} dx"' , (4) 

dS = pady" - Hadx"" , (5) 

where we used Poisson brackets. With equation ([5]) we are able, if desirable, to build an explicit 
solution for the HJPDE, which gives the family of surfaces orthogonal to the characteristic curves. 

3 Integrability 

The vector fields Xa form a basis of vectors in the tangent space of the family of surfaces, as it can 
easily be seen by the HJPDE Xa{S) = XadiS = XaPi = 0- Foi' integrability we need that these vectors 
form a complete set of linear independent (LI) Hamiltonian vector fields. The necessary and sufficient 
condition is given by the Frobenius' integrability condition: the Lie derivative of a member of the base 



2 



with respect to the Hamiltonian flow generated by other member must be a vector in this tangent 
space. In other words, they must close a Lie algebra: 

[X,,X^] = C^/X^ . (6) 

Since the constraints are the generators of the flows, this condition is reflected on the system of HJPDE 
which must form a system in involution: 

{,/.„,<A/3} = C,/0, . (7) 

However, systems in physics do not, in general, respect this condition only with the constraints 
extracted from the Lagrangian. The problem rests in two situations: it may happen that the set of 
HJPDE is not complete. In this case the integrability conditions should provide new constraints to 
complete the system. This scenario is actually well supported by the conditions in the form ([T]). Of 
course, if the final complete set of constraints closes a Lie algebra the system is completely integrable. 

The other situation, which is not covered by condition ([7]), is the case in which the vectors in 
the tangent space of the family are not really LI. To deal with this case we will use as integrability 
conditions the equations 

#a = Ca/3^</'7 = , (8) 

with use of the constraints equations. However, it must be noticed that ([7]) is covered, since d(j)a = 
{(f>a,4'f3} dx^ . This is important since the Frobenius' condition is a theorem and must be satisfied for 
every integrable system. 

The conditions ([8]) can lead to three situations. It can result in expressions of the type f{y, x,p, vr) = 
and these expressions are the ones we are interested in first place, because they must be considered 
as constraints in equality to the former set. They must obey the integrability conditions as well and 
we should be aware of the appearance of other constraints. When all possible constraints are found 
they must be inserted in the formalism. 

However, there will be no independent variables related to the constraints that come from the 
imposition of integrability. We will deal with this problem by expanding the space of independent 
variables, relating arbitrary parameters to each new constraint. The set of all constraints, (pz in which 
the index z now cover all the expanded parameter space, is supposed to be complete and the new 
characteristic equations can be derived from the fundamental differential 

dF = {F, dx" , (9) 

when, now, is the set of all independent variables, including the ones related to the new constraints. 

At this stage we must use the differential ([9]) to test the integrability of all constraints again. It 
can happen that a subset, if not all the constraints, identically satisfies the conditions. This is the 
second situation that will happen when this set is actually in involution, obeying the condition ([7]). If 
all HJPDE are in this case the system becomes completely integrable. If only a subset is in involution, 
the system is only partially integrable. 

4 The Generalized Brackets 

The first two possible results of the integrability conditions ([8]), related above, are in agreement with 
the conditions ([7|). There is, however, a third possibility not predicted by the latter. It may happen 
that the final set of constraints turns out to be not in involution with the PB. 

One way to deal with this problem, using d?]), is expanding the space of constraints finding new 
vector fields that closes, finally, a system in involution [9j. However, the formalism itself cannot give a 
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consistent way to find these constraints: any system that leaves the final set in involution is a good one. 
Other way [TO] is performing a canonical transformation that takes pairs of non-involutive constraints 
into conjugate variables. 

The conditions ^ provides a more natural way to deal with non-involutive constraints. It permits 
linear dependence (LD) of the tangent vector fields X„. 

Let us suppose that all constraints are found and form a complete set of HJPDE, = 0, which are 
related to a set of independent variables x"^. Separating the time variable from the other parameters 
the integrability conditions ([8]) reads 

#a = {4>a, 4>o} dt + {(pa, 4>z} c/x^ = . 

First, WG s-iicilyzG the conditions for (j)x'- 

{(px, (po] dt + {</)^, (pz} dx^ = . 

Let us define an antisymmetric matrix from the PB between the constraints (pz, Mxz = {<Px,4'z}- 
Then we can write 

Mxzdx' = -{(t>x,Mdt . (10) 

If the complete system of HJPDE is not in involution, we deal with a regular M matrix. In this case 
it becomes clear that all vectors Xa are LD from the total differential equations 

dx' = -{M-^Y^ {cPxAo} dt . (11) 

Notice that the fundamental differential ^ becomes 

dF = {F, (Pa] dx"" = {F, ^o} dt + {F, 0J dx^ 
= [{F, (f>o} - {F, (Pz} (M-T W., <Po}] dt . 

Therefore, we redefine the dynamics by eliminating all the independent variables with exception of 
t. It allows us to define the Generalized Brackets 

{F, G}* ^ {F, G} - {F, (Pz} {M-'r {^^, G} , (12) 

which have all the properties of the PB: it is a bilinear antisymmetric operator that obeys the Jacobi 
identity and the Leibniz rule. With the GB the dynamics is given by 

dF = {F,^o}*dt . (13) 

Besides, all the constraints are in involution with the GB, {(j)a, (f'ls}* = 0- There is no need to test the 
integrability of (po, we can easily see that it is identically satisfied with the GB. 

Let us suppose the case in which the M matrix is singular of rank r < k. In this case there will be 
a set of r non-involutive constraints. The equations (jlOh can be written by 



Mxadx^ + Mxzdx' = -{(px, (po} dt , 

where {a} = {1, . . . , r} and {z} = {r + 1, . . . , k}. There are two sets of equations: 

M-^^dx~^ + M-^,dx' = -{(Pi,(Po}dt (14) 
Mx-adx^ + Mx-zdx^ = -{(Px,(po}dt . (15) 

The first set gives 

M-badx^ = - {05, M dt - {(ph, (t>z} dx' = - {^i, (pa} dx^ , 
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making {a} = {0, r + 1, . . . , k}. Since the rank of M^^ is r, Mj^ is regular and we can write 

dx^ = -{M-^f^{c^i,c^a}dx° . (16) 

Hence, we can eliminate r independent variables. The fundamental differential becomes 

dF = {F,(Pa}dx'' + {F,(P^}dx'' 

= {F, dx^ - {F, (M-'f {<!)-„ <j)^} dx^ 

= {F,(^^rdx^, (17) 

in which we have the GB 

{F, G}* ^ {F, G} - {F, <p,} {M-'f {05, G} . (18) 

The set of equations ([15]) gives {(px, (pa}* dx°' = 0. If Mxz has rank r, the structure of the GB gives 
{4'x,4'z} = and {(px,4'a} = 0, which left us with the conditions 

{(l)xAoV = 0. (19) 

In the beginning of this analysis we made the supposition that all the constraints that could come 
from integrability conditions are computed in (f)a. In this case the above conditions are just — 0, 
identically satisfied. But we could define the GB for any set of non-involutive constraints. If there is 
any new constraint yet to be found it will appear in conditions (jl9p . 



5 Free particle on a surface 

In this example we will apply the described procedure in detail in order to clarify the main aspects of 
the HJ formalism. Let us consider a free particle in an Euclidian space M" restricted in a m dimensional 
surface M."^. The surface is defined by the set of equations "02 (x) = 0, in which {z} = {1, . . . , k} and 
k + m = n. We will work in the following notation: x represents a dot in whose coordinates are 
given by the set (x^, x^, . . . , x"). All vectorial quantities are represented in the same way. 
The Lagrangian can be written as 

L(x,v) = i mv2 + n^V;,(x) , (20) 

in which v = x is the velocity of the particle and are Lagrange multipliers. The momenta are given 
by 

dL dL 
P = = "IV , vr^ = — = . 

We have by these relations n velocities and k conditions over the momenta. The canonical Hamil- 
tonian function of this system is given by 

Ho = p-v + TT^u^ -L = ^- u'i)-, . (21) 

Im 

Therefore, we have the following set of constraints: 

00 = vro + = , (22) 
0.;1 = TT, = . (23) 
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We should analyze the integrability of this set. The evolution of any dynamical function F(x, p, u^, tt^) 
is given by 

dF = {F,(Po}dt + {F,(j),,,}du' . 
Because {(pz;i,4'x;i} = 0, the condition for 4>z;i, 

d(t>z;l = {<l>z;l, Mdt + {(pz;l,(px;l}du'' = , 

gives a secondary set of constraints 

4>z;2 = tpz=0 , 

which are expected since those are the equations of the surface. These constraints must also obey 
integrability. Hence, we apply 

#2;2 = {4>z;2, + {(t)z;2, (l>x;l}du'' = , 

which gives a tertiary set: 

(l>z;3 = P • VV-z = . 

These constraints tells us that the momentum p is tangent to the surface. And, from these, 

d(t>z;3 = {<l}z;3, (t>o}dt + {(t)z-3, (t>x;l}du'' = , 

a quaternary set arise: 

(l)z;i = mu'^A^z + p • V (p • Vi^z) = , 

in which we define the matrix Axz = Vi/'o; " VV'z- 

Notice that {(j)z;4:, (po} 7^ and {4>z;i, 4'x;i} 7^ 0, so that the integrability of 02;4 will not give any new 
constraint, but a total differential equation relating the independent variables. We suppose, therefore, 
that all constraints are found and the complete set is given by 

(t>Q=Po + Ho , 

4>z;l = T^z , 

4'z;2 = i^z , 

4'z;3 = P • "^Ipz , 

^2;4 = mAzxU'' + p • V (p • Vlpz) ■ 

To consider this set in the theory we must build a differential that contains these constraints as 
generators of the dynamics. However, we do not have independent variables related to all constraints in 
the system. Therefore, we will expand the parameter space with new arbitrary independent variables 
{uj^, T^,0^). Then we define the dynamics of the system by the new fundamental differential 

dF = {F, M dt + {F, ^z;i} du' + {F, </>,;2} dio' + {F, (j^z-s} dr' + {F, <t>z;A} dO' , 

as prescribed in section 3. 

The above set of HJPDE is not in involution, as we can see with the PB 

{(f^x-A, (t>y;2} = {4>x-l, (t>y;3} = 

{(px;l, (t>yA} = -mAxy {^a;;2, 0y;3} = ^xy 

{(t>x;2, (f>y;4} = P • V {A^y) {(t)x;3, 0j/;4} = ^xy 

where 

Txy = -mu'V {Ayz) • VVa; + P • V [p • V (A^^)] - 3 (Wa;) • [p • V [p • V {V^Py)]] . 
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We expect, therefore, that all independent variables can be eliminated in function of t. Let us reorganize 



the constraints so that their order in the M matrix becomes 
by (/.J2, = {1,4,2,3}) 



')z;l,(Pz;4,(Pz;2,(Pz;3) 
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The matrix is given 



(24) 



Let us suppose that the surface is smooth, of class C°°. This condition guarantees the existence of 
the inverse 



(M" 



where 
and 



/ 1/m {A-'^yy T'^y 

-l/m (A-^)^'s^ 

_fxy Q Q 



^xy ^ (A-i)^'^r^^(A~i)"'2/ 



j^xy 



-(A-i)^s^ 




(25) 



j^xy = p . v(A 



-l\xy 



This inverse is unique. Since the matrix is regular all independent variables are indeed dependent of 
the parameter t, as learned with equation (llip . We are able to define, then, the GB of the system: 

{F,Gr^{F,G}-{F,<t>jJ{M-'y^'y{4>j^,G} . (26) 

The only non zero fundamental GB are given by 

{^,pr = I-VM^~T^i^y^V . 

All others are zero, including {u, vr}*, which is expected since are degenerate variables in the theory. 
Here we define the projector operator V, whose job is to take vectors in R" into vectors on the surface. 
It is a singular operator, whose base of the null space are the vectors Vipz- This is expected since 
normal vectors are projected into null vectors on the surface. 

Let us go for the characteristic equations of the reduced system: 



dx = {x, (/>o}*dt = — V ■ p dt 
m 

P = V-p, 



(27) 



dp = {p, (poYdt = V ■ u'^Vtpzdt = 
a = V ■a = . (28) 

The first equation gives the projection of the momentum vector over the surface, it is a geometric 
relation. The second equation is the dynamical equation: it tells us that the acceleration induced on 
the surface is zero. Of course, we cannot invert the projector, since it is singular, so the equations of 
motion are still degenerate. Since p = 0, the derivative of (I27p over time is given by 

^=^.i. (29) 
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Let us suppose the case in which n = 3 and the surface is just with equation ^ = — = 0. 
We have VV' = 2x and A = 4x^. In this case A~-^ = l/4r^ and 

where (8> represents a dyadic product. Using the fact that (pz-.s = imphes x • x = 0, equation (|29p 
gives 

i + (j^x = , (30) 

in which uj = v/r and v = |v|. The case in which r = 1 and m = 1 is an analogous to the hnear 
quantum mechanical u-model treated in |10j . 



6 The multi-dimensional rotator 

Let us work with the Lagrangian |10j 



L = ^ v2 + MX • V , (31) 

where v = x. This system constitutes a more general case than the previous example. The conjugate 
momenta are given by 

p = V + ux , 
TT = . 



The canonical Hamiltonian gives 



Ho = ^[p- u^f = ^ P' , (32) 



in which we define p = p — ux. Therefore, the system has the constraints 



= Po + Ho, (33) 
H = TT . (34) 



The fundamental differential is, in this case, 

dF = {F, (l)o}dt + {F, (l)i}du . 
The integrability condition for c^i gives 

d<j)i = {(pi, (j)o}dt + {(pi, (pi}du = p • xdt , 



which results in a new constraint 



^2 



p • X = p • X — nx^ = . (35) 



The integrability of (p2 gives no new constraint. The system is not in involution, giving a regular M 
matrix and its inverse 

- ' -x2 j ' M ^^^2 P 
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The GB of this problem, that will eliminate the independent variables in function of t, is given by 
{F, G}* = {F, G} - {F, 4>i} (M-i)i2 {^,, G} - {F, <p,} {M'^f^ G} , 
which gives the following GB: 

{x,p} = I 2~ = ^' 

{y.,uy = ^ , {u,p}* = ^ . 

x^ 

The characteristics equations are given by 

dx = {x, (j)Q}*dt = V ■ p dt , 

dp = {p, (/)o}*dt = , 
p2 

du = {u,(l)o}*dt = dt . 

x^ 

Hence, we have the equations 

V = V-p, (36) 
^ = 0. (37) 

Derivation of the first equation over time and the use of the second equation gives 

V • p 



± = V-p 



X2 



where we used the constraint p • x = 0, which also implies, from (|36p . v = p. We have, therefore, the 
nonlinear second order ODE 

i • X + = . (38) 

This equation can also be written by 

*^ = °' 

Notice that the equation of motion (j38p is obeyed by any system in which x • x is a constant. The 
most simple particular solution is given by |x| = ^/at + f3. However, considering a system in which 
x^ = is fixed and |v| is constant, we can build a solution 

X = a.smuj{t) + h cos uj{t) 

where a and b are constant vectors and 

co{t) = ^{\^r\t + f3), 

where /3 is a parameter. This solution is of the same class of the one found in [TO] for the n dimensional 
case. 
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7 Landau model 

Consider the Lagrangian [10 



L = ^ {mviVi + BeijXiVj - kxiXi) {%} = {1, 2} , 

which represents a charged particle in a plane with a transversal magnetic field B and a harmonic 
potential. In the limit m = this system is described by 

L = ^ {BeijXiVj - kxiXi) . (39) 



The conjugate momenta give the relations 



Pi — 2 ^^^i^i 



Let us write the canonical Hamiltonian: 



-f^o = ^ kxiXi , (40) 



that gives us the system of HJPDE 



IJo + ^0 = , (41) 
1 
2 



Pi + - BeijXj = . (42) 



This problem is completely constrained and all x variables are parameters in the theory. The 
fundamental differential is 

dF = {F,(pQ}dt + {F,^i]dx' . 

Since {(j)i, (pj} = Bcij, the constraints are not involutive and the integrability conditions for (pi will not 
give any constraint, but a total differential equation relating the variables x and t. We have the M 
matrix and the inverse 

M,,=Beij ^^^ = --^ ■ 

The GB are defined by 

{F, GY ^ {F, G} + i {F, 0j6,,{</)„ G} . (43) 

Therefore, 

{Xi,Xj"\ = — Cij , {^i^Pj} ~ 2 ^ij ' \PiiPj"\ ~ ^ ^^ij • 

The equations of motion are given by 



dxi = {xi,(j)o}*dt = —— keijXjdt , 



1 

B 
1 



dpi = {pi,(t>o}*dt = -- kxidt . 
The result is 

^2 

+ Xi = , (44) 

in full agreement with [lOj . 
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8 The Giiler's example 

Consider the following Lagrangian function ^ 

^ = ^ - ^ {(12- qsf + bq2-c , (45) 
in which b and c are functions of qi, q2, qs and t. The momenta are given by 

Pi = qi, P2 = -^{q2- qs) +b , P3 = ^{q2-q3) ■ 
Therefore, there is a constraint relating the momenta p2 and p^: 

P2+P3-b = . 

We will choose to write the Hamiltonian in terms of pi and p2, which implies to choose qi and q2 as 
dependent variables, while qs ant t are the independent ones. We have, then, 

Ho = ^pj-{b-p2f + c, (46) 

which gives us the set 

(1)0 = Po + Ho, (47) 

(Ai = P2+P3-b = . (48) 

The fundamental differential is given by 

dF = {F,(j)o}dt + {F,(t>i}dq3 . 

With this differential we test the integrability of 0i, which gives a new constraint 

„ , db db db dc dc 

dqs dt dqi dq2 dq3 

With b = qi + qs and c = 51 + 52 + ^3 this constraint can be written as 

(p2^2p2+pi-2qi + l = . (49) 

We have {(j)2,(l>i} = 1 and no new constraints. The M matrix and its inverse follows: 

/ -1 \ / 1 

^=(1 J ' ^ =(-1 

The GB are defined by 

{F, G}* ^ {F, G} - {F, <Ai}{02, G} + {F, MWuG} (50) 

and gives the following non zero fundamental GB: 

{^1,92}* = {^1,93}* = {gi,P3}* = -1 

{92,93}* = {92,^3}* = {pi,P3}* = -2 
{92,^2}* = {q3,P3}* = 1 
te,Pir = 2 . 
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Those give us the equations 



dqi = {qi,(t>o}*dt = {2qi - 2p2 - 1) dt 
dq2 = {q2, M*dt = (4gi - 4p2 + 1) dt 
dqs = {q3,M*dt = {2pi+2p2-2qi-2q3 + 3)dt 

dpi = {pi,(po}*dt= {Aqi+2q3-4p2)dt 

dp2 = {P2,(t>o}*dt = -Idt 

dp3 = {P3, M*dt = {2pi - 2q3 + 3) dt , 

which can be written, using the constraints, by 

qi = 2 (gi + 1) , q2 = 4 (gi + 1) , Q3 = 4 (^(73 - ^ • (51) 

As we can see, the solutions in [9] satisfies these equations. 



9 Final Comments 

The appearance of non-involutive constraints violates the Frobenius' integrability condition, by which 
the constraints related to a set of LI Hamiltonian vector fields must close a Lie algebra with the 
Poisson brackets. In this work we have shown that the imposition of the integrability conditions 
d(f)a = provides a self consistent way to treat non-involutive (second class) constrained systems. 

It must be stressed, however, that the procedure do not change the Frobenius' theorem: if a system 
is in agreement with ([8]) it is also in agreement with ([7]). We saw that the imposition of ([8]) can reveal 
new constraints or show that the system, or part of it, is integrable, cases that are covered by ([7]), 
but it can also reveal linear dependence in the vector fields tangent to the surfaces, which implies 
elimination of independent variables related to non-involutive constraints. 

The process of elimination shows that it is always possible to redefine the dynamics with the 
introduction of Generalized brackets. The GB actually defines a symplectic structure in a reduced 
phase space. Systems that presents only non-involutive constraints (M regular case) are completely 
reduced and the evolution of any observable is generated by the canonical Hamiltonian with the GB, 
as shows equation (|13p . Therefore, non-involutive constraints are not dynamical generators in the 
reduced phase space. This is reflex of the fact that the constraints that were not in involution with 
the PB are in involution with the GB now. 

However, if the system presents secondary constraints a previous step is needed. Since there 
are no independent variables related to constraints that come from the integrability conditions, we 
may expand the parameter space to include new parameters, one for each of those constraints. The 
secondary constraints are, then, assumed to be generators of infinitesimal transformation^ in the phase 
space, in equality with the primary constraints, as shown in the evolution equation ([9]). 

If the system has only a subset of non-involutive constraints we have the case in which M is singular. 
This subset of r constraints is computed in the GB of the system. However, integrability conditions for 
the other k — r constraints are given by equations (jl9p . which can only result in new relations between 
the variables (in this case new constraints), or in identical = relations. The primary case requires 

^This statement can always be made, but we stress that secondary constraints may not be dynamical generators in 
the reduced phase space. This is actually the case with non-involutive constraints. For involutive set of constraint the 
validity and consequences of this statement are still under investigation. 
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the insertion of these new constraints in the system. However, this case wih not occur if all secondary 
constraints are already previously found. The second case implies that the remaining constraints must 
be in involution with the GB and the system whose evolution differential is given by (llTh is already 
completely integrable. 
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